NAVIER STOKES EQUATIONS ON THE /3-PLANE 



M.A.H. AL-JABOORI AND D. WIROSOETISNO 

Abstract. We show that, given a sufficiently regular forcing, the solution 
of the two-dimensional Navier— Stokes equations on the periodic /3-plane (i.e. 
with the Coriolis force varying as /o + 0y) will become nearly zonal: with 
the vorticity u)(x,y,t) = ui(y,t) + u>(x,y,t), one has |w|h- s < P~ 1 M S (---) as 
t — > oo. We use this show that, for sufficiently large /3, the global attractor of 
this system reduces to a point. 



1. Introduction 

The two-dimensional Navier-Stokes equations (2d NSE) have been the subject 
of many studies and its basic mathematical properties (existence, uniqueness, reg- 
ularity, etc.) are now well understood; see, e.g., [5] [11] for reviews. As a tool to 
understand various geophysical flows, it is often desirable to include the effect of 
planetary rotation, but a constant rotation rate (the so-called /-plane approxima- 
tion) has no effect on the dynamics when periodic boundary conditions are used. 
To feel the effect of rotation, we need to go to the so-called /3-plane approximation, 
in which the rotation is given by fo + /3y. 

Simple physical arguments and numerical studies [15] suggest that a rotation 
rate that varies as j3y tends to force the solution to become more zonal (a zonal 
flow is one that does not depend on x). In this article, we prove that this is indeed 
the case, by obtaining a bound |w(£)|^ s < (3~ 1 M s (f; • • • ), valid for large time i, on 
the non-zonal part u of the flow in terms of the forcing /. 

With the further assumption that the forcing is independent of time, it has 
been shown that the Navier-Stokes equations possess a global attractor A of finite 
Hausdorff dimension. The long-known and nearly optimal bound on this dimension 
[4] also applies to our rotating case, but it does not take into account the effect 
of the rotation. Using our bounds on w, we show that the dimension of A is zero 
for e sufficiently small, reducing the long-time dynamics to a single steady (and 
stable) flow determined completely by the forcing. This is to be contrasted with 
the situation for larger (but still small) e, where the solution, although nearly zonal, 
evolves in time even though dtf = 0. 

Among the works similar in spirit to the present article, we mention [5] where 
weak convergence to zonal flow is proved for the (more difficult) /3-plane shallow- 
water equations. A related result for the inviscid Euler equation can be found in 
[lOj . The technique of using rapid oscillations to obtain better bounds have been 
used in different contexts in, e.g., [21 E]. 
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In the rest of this section, we describe the problem and set up the notation. 
In Section [21 we review basic results on 2d NSE which will be needed later. The 
heart of this article is Section [3J where L 2 and H s bounds are obtained for the 
non-zonal component of the flow. An application of these bounds to the dimension 
of the global attractor follows in Section |4j The proof of an L°° Agmon inequality 
is presented in the Appendix. 

In dimensional form, the two-dimensional Navier-Stokes equations read 

(1.1) d t v + v ■ Vv + fiyv 1 - + Vp = p,Av + /„ 

where the constant rotation fo has been dropped since it has no effect (i.e. in 2d 
NSE, there is no difference between equatorial and mid- latitude /3-planes). Here 
v = (u,v) is the velocity with v := (— w, u) and p is the pressure obtained by 
enforcing the incompressibility constraint V-w = 0. In what follows, we will work 
with the dimensionless form 
Y 

ft o\ d t v + v-Vv-\ v 1 - + Vp = p,Av + f , 

(1.2) e 

Vv = 0. 

We work with x = (x, y) G M, := [0, L\\ x [— i 2 /2, ia/2], with periodic boundary 
conditions in both directions. Note that we have replaced f3y in (|1.1[) by Y(y)/e, 
where Y(— L2/I) = L2/2 and Y(y) — y for y e (— L2/2, L2/2]. Furthermore, we 
assume the following symmetry on the velocity 

(1.3) u(x,-y,t) = u(x,y,t) and v(x, -y, t) = -v(x, y, t). 

It is readily verified that if the initial data v(0) and the forcing f v (t) also satisfy 
this symmetry, which we henceforth assume, it persists for all t > 0. Note also that 
periodicity and (|1.3[) imply that 

(1.4) v(x, -L2/2, t) = v(x, La/2, t) = 0. 

With no loss of generality, we require that the integral over M of v vanishes. 

In two dimensions, it is convenient to work with the vorticity uj := V -v = 
d x v — d y u, whose evolution equation is 

Y' 

(1.5) d t u) + v ■ Vw + — v = p,Au + f. 

Here / := V ■ f v and we can recover the velocity using v = V^A -1 ^. By our 
assumption on v, the integral of u> over M is zero; similarly, A -1 is defined uniquely 
by the zero-integral condition. The symmetry (|1.3j) implies that u){x,—y,t) = 
—uj(x, y, t) and 

(1.6) u(x, -L 2 /2, t) = u(x, L 2 /2, t) = 0. 

Now Y'(y) = 1 — L2 S(y — L2/2), where 6 is the Dirac distribution. Using the fact 
that v(x,±L 2 /2,t) = 0, we replace vY' by v in (|1.5p and write 

(1.7) d t uj + v -Vcj + -v = fiAui + f. 

This is the form that we will be mostly working with. 

It is also convenient to write (ll.7[) in the usual functional form 

(1.8) d t u + B(u,u) + -Lu + pAu = f, 
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where B(uj,J) := (V^A"^) • Vw», Luj := dxA^u and A := -A. The follow- 
ing properties, valid whenever the expressions make sense, are readily verified by 
integration by parts and the boundary conditions (|1.4|l - (jl.6|l 

(B(w,J),J) L2 =0, 

(1.9) (Loj,oj) L 2 =0, 

(Auj,uj) L 2 = |Vw|£2. 

2. Preliminary Estimates 

The estimates derived in this section are standard from the theory of 2d NSE 
(see, e.g., [121 El IB]), with very minor modifications to handle the Coriolis term. We 
gather them here for later use. 

We start by noting that the vanishing of spatial integrals of v and to implies the 
equivalence of the norms \uj\ hs and |V s w| i2 := |(— A) s / 2 w| i2 , which will thus be 
used interchangeably below. We denote by Co the constant in Poincare inequality 

(2.1) c \V s ujy < \V s+1 u\ L2 . 

Besides the usual Sobolev and interpolation inequalities for two and one dimensions 
(for functions depending on y only) , we note the one-dimensional Agmon inequality 

(2.2) H^cl^lVu,^ 2 . 

A version we use for the two-dimensional case is in Appendix |XJ 

The L 2 estimate for the velocity is obtained by multiplying (|1.2|) by v and using 
Cauchy-Schwarz, 

(2-3) ^\v\h+n\Vv\l 2 < C -\fv\l>- 

Assuming that /„ G L^L 2 ,, we thus have v G L^L 2 , CiL 2 !!^.. Here and henceforth, 
L^L^, := L p ((0, oo); L q (M)), and H x and |V s w| below are defined in the usual way. 
We denote 

(2.4) [to] := sup \w(t)\ L 2. 

t>o 

Here and elsewhere in this article, c denotes a generic constant depending only on M 
whose value may not be the same each time it appears, while numbered constants 
such as c Y have fixed values. 

Now let (p G L°° be such that ip'(t) = 1 for t G [0, \] and tp(t) = 1 for t > 2, so 
ip(t) ~ tanhi. Multiplying (| 1 . 7[) by ipuj in L 2 , or equivalently, multiplying (|1.2b ) by 
—ipAv in L 2 , we obtain 

(2.5) ^(pMiO+M^lVwIi, <^'|w| 2 + ^|/Ji 2 . 

Assuming henceforth that [/„] < oo, we have uj G L^L 2 , n ^ 2 -ff^. and, for t > 

t (Ko)| l2 ,[/j ;m ), 

(2.6) \u(t)\ L3 <-I/J. 

M 

Note that To does not depend on u>(0) and that the requirement f v G LfL 2 x can 
be weakened in t, but we shall not do so here. 
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A bound in H m is obtained as follows. Fix a multi-index a — (0:1,0:2) with 
\a\ := ai+a 2 = m, and multiply (dHJ) by D 2oi lo := d% ai d% a2 u) in L 2 , 

(2.7) ~\D a u\ 2 + (v • Vw, D 2a Lu) + -(«, D 2a oj) + /1 |VD Q w| 2 = (/, D 2q w) 
2 dt e 

where here and henceforth | • | and (•, •) denote L 2 norm and inner product. The 

linear term involving 1/e vanishes, and one then proceeds as usual: Using the fact 

that (v ■ VD Q w, D°w) = 0, the nonlinear term is bounded as 

I (v ■ Vw, D 2a u) I < I (( D/3?; ) ' VD a ^w, D°w) I 

1<I/3|<|«| 

(2.8) < ID^I^ID^+^I^ID^I^ 

< c(m) E£i |w| H i_ 1/a |o;| Hm _ 1+ 3/ 3 |w| Hro 

where we have used Sobolev inequalities for the second and third line, and where 
I := |/3 1 in the last line. Using the interpolation inequalities 

I , I ,(2m-2/+3)/(2m+2) | , (2/-l)/(2m+2) 

^ " ' L, I ^ .1, ,|(2i-l)/(2m+2), |(2m-2*+3)/(2m+2) 

followed by Cauchy-Schwarz and summing over a, we obtain 

(2.10) ^\ 2 Hm + f \^ Hm+1 < ^ \u,\Uu\ 2 Hm + ^ \f\% m ^ 

for to = 1, 2, • • • . Proceeding by Gronwall and induction on (12.10[) . we have the 
following uniform bounds independent of the initial data 

c(TO) . „ N m + l 



Hr)\ 2 Hm+1 d r<*(^lV m -'ff + l) m+1 
fj, \ n / 

valid for all t > T m (\v(0)\ L2 , [V™ 1-1 /] ; /i). Here v — c 2 fi with c the constant in 
Poincare inequality (|2.1j) . Note that for T m to depend only on |v(0)| i2 , and for the 
validity of (|2.11[) for v(0) ^ H 1 , we need to multiply by ip as in f|2 .5(1 . but this was 
not done explicitly for conciseness. 

3. Bounds on the Non-zonal Component 

Assuming sufficient regularity for /, which implies that for u for any t > 0, we 
expand them in Fourier series 

^*) = E fc ^)e ifc -*- i0fct/£ 

/0M) = E fe /fc(*)e ife '* 

where fc = {ki,k 2 ) G Z L := {(27rfi/Ia, 2nh/L 2 ) : (h,l 2 ) <= Z 2 } and Q k ■= -h/\k\ 2 
is i times the eigenvalue of the linear operator L for wavenumber fc. Since uj and 
/ have vanishing integrals over M, ujq = and /o = 0. Here and in what follows, 
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sums over wavenumbers are understood to be taken over Zl. In terms of Fourier 
components, (|1.7|) reads 



(3-2) ^ + E jk B juUj u h e^-^-^We + ^ ^ = fi ^t/e 

where the coefficient of the nonlinear term is 

(3.3) B jkl = (BCeW-.e 1 *-).^-) = |M| ^*i+fc-i 

with j A k := jifc2 — J2&1. We note that the linear term e^Llo has been removed 
from (|3.2j) by including exp(— ififei/e) in p. lb ). 

Let us split to into a slow part a), for which fl^ = 0, and the remaining fast part 
Co := to — Co, viz., 

*(*,*) =E fcl =o^(*)e ife - a! 

We note that, also having zero integrals over M, tl> and Co are orthogonal in if m for 
m = 0, 1, • • ■ . For convenience, we also define 

/n rN _ Jwfc if fci = _ Jwfc if ki ^ 

(3.5) Wfc := < and lo^ := < 

I otherwise, I otherwise. 

Our objective in this section is to obtain long-time bounds for Co that tend to zero 
as e — > 0. 

3.1. Bound in L 2 . The development in this subsection largely follows that in [13] 
for the primitive equations, the main difference being the absence of a spectral gap 
(that is, the eigenvalues of the antisymmetric operator L accumulate at zero in the 
present case). 

We start by multiplying (11.81) by Co in L 2 , 

(3.6) (8tw, u) + (B{lo, u), Co) + ^(Llo, Co) + n (Aw, Co) = (/, Co). 

Now, using (|1.9b ) twice and the fact that B(Co,Co) = 0, 
(B(lo,lo),Co) = (B(lo,Co),Co) + (B(lo,uj),Co) 

(3.7) = (B (Co, Co), Co) + (B (Co, Co), Co) 
= -(B (Co, Co), Co). 

Thus p.6[) becomes 

(3-8) ~\Co\ 2 + fj, \\7Co\ 2 = (B(Co, Co),Co) + (/, Co). 

Dropping the nonlinear term for the moment, the fact that Co is rapidly varying 
while / is slow implies that the effective forcing from the rhs becomes weaker for 
smaller e. This essentially is the mechanism for the attenuation of the fast part 
Co; the nonlinear term will be handled in the proof below. Recalling the definition 
(|2.4[) . we state the result of this subsection. 
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Theorem 1. Assume that the initial data v(0) £ L 2 (M) and that the forcing is 
bounded as [V 2 /] + fdtfi < oo. Then there exist To{\v(0)\ L2 , [V 2 /], fi) and 

M ( [ V 2 /] , \d t /] ; ^ ) «ucA that, for t >%, 

m)\ 2 L 2<eM , 

(3.9) , /•«+*' 

/ie-"< f+f ' / e^|Vw(T)| 2 dr<eM . 



Proof. Recalling that v = /ic 2 ,, we obtain from ([3 

(3.10) ^-ie vt \Cj\ 2 ) + ^ ut \VCj\ 2 < 2e^ t (B(w,£),cj) + 2e l/ *(/,w). 
at 

We integrate the last term from to t by parts, 

[ e» T (f,Q)dr = |M|Efc / /feW^T)e i!2 ^/ £ +" T dr 
Jo Jo 

where the prime on the sums indicates that the resonant terms (i.e. those with 
flk = 0) are excluded. Defining the operator by, for any w for which it makes 
sense, 

d* t w :=e~ tL / £ d t (e tL/E w) 

(3-12) 1 

=>■ dfib := dtuj H — LGj = —B(oj, d) — /iAu> + f, 

£ 

and defining the operator In by 

(3.13) !„/>, t) := £ k ^ /fc (t) e- = i £ ^ /*>(*) e ife *, 

which being the restricted inverse of L is also antisymmetric, we can write 

f e^(/,w) dr = e(l n /,w)(t)e ,rt -e(ln/,w)(0) 

(3.14) - 70 t 

-e / [Klf 2 /,^) + (9 r l f2 />) + (lo/,9^)]e^dr. 
Jo 

Using (|3.13[) . the endpoint terms can be bounded as 

(3.15) |(ln/,0)|<c|V/||Vw|. 

We now bound the terms in the integrand. First, 

(3.16) \{dr\ a f,S>)\ = \{d T f,\ n u)\ <c|0 T /||Aw|. 

Next, using (|3. 12b ) and noting the fact that (In/, /) = 0, we bound the last term 
in flQlH by 

(3.17) |(ln/~MAw)| <^c|A/||Aw|; 
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and, using Sobolev and interpolation inequalities, 

\(\ Q f,B(oj,uj))\ <c\Vf\ L2 \VB(Q,Q)\ L2 
(3.18) < clV/liiM^lVwI^ + clV/l^lV-^I^IAwl^ 

<c|V/||Vw||Aw|. 
Thus the integral in (13. 14)) is bounded as 

' Kb/» + (9 T b/» + (In/, TO] dr 



(3.19) 



< c /" [/x | A/I |Aw| + |fl T /| |Aw| + |V/| |Vw| \Auj\]e VT dr 
Jo 

<c[ {(1 + M ) |A£| 2 + /u | A/| 2 + |9 T /| 2 + |Aw| |Vw| |V/|} eT dr. 
Jo 



We now treat the penultimate term in A3. 10[) . First, we write 

(B(w,w),w) = E ifeJ B jM u j u k uie- l ( a * +n *W e 
(3.20) i 

and then note that Bjki + B^ji = in the resonant case, i.e. when Qj + il^ = 
and l\ = 0. Furthermore, we have 

(321) ^ + ^-(w + w)™ = {iAk) iw-w) m 

= jilj(ijja-4j) M =-I 2 (!! J + !!»)|M| 

whenever j + k = I and Zi = 0. Motivated by (|3. 201) . we introduce the bilinear 
symmetric operator i?o by 



(3.22) 



* — ^-ifcZ J 



2i ^jfcJ 

for any w", and w, where the prime on the sum again indicates that resonant 
terms (for which fL + fi^ = 0) are omitted. We note that, thanks to (|3.21l) . the 
resonant terms are also absent in (B(ui,ui),ui). Integrating by parts, we have 

e VT (B(w,w),u>) dr = ee w (Bn(w,tj),u))(i) - e(B n (w,w),w)(0) 



(3.23) 



c VT [ v {B n {u,u),oj) + 2(B n (.d* T u,u),uj) 



+ (Bo(a),d)),a T w)] dr. 

For the last term in the integrand, we use the fact that B(Cj,uj) — B(u>,Q) 
B(ui, ui) = to write 

(3.24) d T ui = -B{Cj, Q) - fiAiD + f 
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and estimate, using H 1 C L°° for / and (|A.2|) for the L°° estimates, 

\(B a (ui,ui),d T u))\ < c\ti\ L2 \d y u\ L2 \f\ L oo + [xc\ui\ Li \d y uj\ L 4 | Aw| i2 
( 3 25) + c |fyw| £a |V -1 w| £ o. | Vw| L 2 

< c |u| |VcD| |/'| + M 1 ' 2 ^ |Aw| 3 / 2 + c |Vw| 3 |w| (log + c ') . 

For the term involving (BQ(d*Li,Li),ui), we bound, using 8*lu + -B(w, w) + /iAw = / 
and the inequality |w| i<x > < c |u> | 1//2 |cZ>'| - 1 / 2 , 

\(Bn(d*u),Q),u>)\ < c\d y f\ \u>\ \uj\ l ^ +c|/| \d y u\ \uj\ l ^ 

+ \xc I Aw I |w| |w'| £ oo + c | V _1 w| L oo | Vw| |w| |w'| L oo 

( 3 - 26 ) < c |V/| |w| |w'| + c l/l |Vw| |w'| + fie p| |w'| 1/2 |Aw| 3 / 2 

-c| w | 2 |Vw| 3 / 2 |w"| 1/2 (log^4 + l V/2 



c |w| 

where all unadorned norms are L 2 . Finally, we bound 
\(B u (ui,ui),uj)\ < c\oj\ \d y u>\ |w| L oo 

< c\u\ |Vw| IcjI 1 / 2 ^'! 1 / 2 . 
Using these also to bound the endpoint terms, the integral in (|3.23[) is bounded as 



(3.27) 



t 



e UT (B(uj,Q),Q)dT 

< ec [\Co\ |Vw| |w| 1/2 |w'| 1/2 ] (i) e" 4 + ec [|w| |Vw| |w| 1/2 |u/| 1/2 ] (0) 
(3.28) r t f 

c | V/| |w| |Vw| + c l/l |Vw| 2 + |w| 1/2 |Vw| | Aw | 3/2 



+ c|w||Va;| 5 / 2 |Aw|^(logi^i [ + c ')} 1 



Putting together ([535]) . (|3~T§)) and (j3T2"gl) . we have 



|w(t)| 2 + p / |Vw| V< T -*> dr < 
Jo 



e- !/t |w( 



+ ec 2 (l+e- 1 '*) sup {|V/||Vw| + |w| 3 / 2 |Vw| 3 / 2 } 

0<t'<t 

(3-29) t 

+ ec 3 (fx) / { |A/| 2 + |9 T /| 2 + |V/| |Vw| |Aw| + |Aw| 2 (l + |Vw|) 



;| |Vw| 5 / 2 |Aw|V^l 0g + /) } e"<-«> dr. 



We now shift the origin of time such that t = corresponds to in (|2.1ip . The 
hypothesis that [V 2 /] + [9*/] < oo then implies that both the endpoints and the 
integral in (|3.29p are bounded uniformly for alii > 0, independently of the initial 
data provided that v £ L 2 initially. Rewriting the bound in (|3.29j) as 

(3.30) \Cj(t)\ 2 + l*f |Vw| V^) dr < e"^|w(0)| 2 + | M ([V 2 /], p t /]; M ), 
the proof is complete. 
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We also note from (f3~29]) that the hypothesis / £ LfHl and St/ G LfL 2 , can 
be weakened to / G L\H* n LfH^ and St/ G I? t L%. 

3.2. Bounds in With a little extra work, H a bounds for w that scales as yfe 
can also be obtained. We do this explicitly for |Vw| and sketch the computation 
for s = 2, 3, • • • . 

For the H 1 bound, we multiply (11.81) by Aw in L 2 to get 

(3.31) ^|Vw| 2 + fi |Aw| 2 + (B(w,w),Aw) = (/, Aw), 
which implies [cf. (|3. 10[) ] 

(3.32) ^-(e^lV^I 2 ) + ijLe vt \ ACj\ 2 < 2e vt (B{uj, w), Aw) - 2e vt (f, Aw). 
As in the L 2 case, we integrate from to t, 



e"*|Vw(t)| 2 - |Vw(0)| 2 + fi / |Aw| 2 e !/T dr 
(3.33) 70 1 

<2/ {(i?(w,w),Aw)-(/,Aw)}e^dT. 
Jo 

The forcing term gives 



(3.34) 



/ e^(/, Aw) dr = e (l n /~, Aw)(t)e I " - e (l n /, Aw)(0) 
Jo 

+ e / [Kln/,Aw) + (S T lo/,Aw) + (l a /,AS;w)]e^dr, 
Jo 

which can be bounded as in the L 2 case, giving 

2 / (/, Aw) e" r dr < ec [|A/| |Aw|] (i) e vt + ec [|A/| |Aw|] (0) 
Jo 



(3.35) 

+ £C / {(l + ^lV^p+^lVVT + lV^/P + lV^UVwIIA/De-dr. 



For the nonlinear term, we use the fact that -B(w, w) = to write 

(3.36) (5(w,w), Aw) = (B(ui,w),Aw) + (B(ui,ui),Aui) + (B (w,w), Aw), 

and, using (£?(w*, w), Aw) = (B(Vw',w),Vu), previously used in (|2.8[> . we bound 
|(B(w,w), Aw)| = |(B(Vw, w), Vw)| < c|w| i2 |Vw|| 4 < c|w| |Vw| |Aw| 

< J|Aw| 2 + ^|w| 2 |Vw| 2 

(3.37) | (B(w, w), Aw) | = |(B(w',w),S y w)| < c |w| L=0 |Vw| 2 2 
|(B(w,w), Aw) | < |Aw| L2 |V _1 w| L oo|w'| L2 < |Aw| i2 |Vw| i2 |w'| i2 

< J|Aw| 2 + ^|Vw| 2 |w'| 2 . 
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Using Poincare inequality on the last term in Q3.37b ). we obtain 

2 [ {B{uj,uj),^Co)e VT At 
(3.38) J ° t 

< c 2 ( M ) J { [loLoc + |V^| 2 ] |V£| 2 + £ |A£| 2 } dr. 

After moving the |Aw| 2 to the left-hand side, a factor of e can be obtained by 
pulling the square bracket outside the integral and using (|3.9b ). Collecting, we 
have 

|Vw(i)| 2 + ^ / t e l/ ( r -*'|Aw| 2 dr<e~" t |Vw(0)| 2 + cesup |A/(t')| |Aw(t')| 

^ 3 - 39 ) +ce f{{\ + m)(|V 3 ^| 2 + |V 3 /T) + |Va r /| 2 + |V^| 2 |A/| 2 } e^-*) dr 
Jo 



ec 



(^Mq sup{|w(OlL» + |Vw(f)r} 



t'>0 

Arguing as in the L 2 case, / G Lf^H^. and <9f/ G LfH^ gives us an 0(y / e) bound 
for ui(t) in L^H^, uniform for large t. 

Bounds in H s can now be obtained inductively. Assuming that Theorem[5]below 
holds for s— 1 (we just showed that it does for s = 2), we multiply (|1.8p by A s Cj and 
integrate the resulting equation in time as above. We bound the nonlinear term 
(B(uj, w), A s Cj) = (B(oj, u), A s u) + {B(lj, Q), A s Co) + (B(u, u), A s u) as follows. The 
first term is bounded exactly as in 



(3.40) \(B(u,u), A s Cj)\ < ^ \V s+1 £o\ 2 + ^ |w| 2 |V s w| 2 . 

1 4 n 

We bound the next term by [cf. (12.81) ]. with |a| = s and 1 < \f}\ = r < s, 
\(B(Q,Q),A S Q)\ < c£ Q/3 |D^| ioc |D Q -^| i2 |D^| i2 

(3.41) <c( S )E: = i|V^||V s -'' +1 ^||V^| 

< c(s) |V s w| |V s w| 2 . 

Finally, we bound the last term as, where now < \/3\ = r < s = \a\, 
\(B(u,u,),A s Q)\ < c£ a/3 |D^| i4 |D Q ^V^| L2 |D^| L4 

(3 - 42) < | |V S+1 ^| 2 + c(s, jtx) |V r - 1 «| 2 / 3 |V r w| 2 / s |V'w| a / 3 |V*- r+1 w| 4 / 3 
< £ |v s+1 w| 2 + c(s, M ) |V s £| 2 |V s+1 ^| 4 / 3 . 

Moving the | V s+1 w| 2 in (|3.40l) and Q3.42p to the left-hand side of the main inequality, 
the right-hand side depends at most on |V S (I>| 2 , which is of 0(e) in L\ from step 
s — 1, and on |V s+1 w| 2 . As before, the worst term (i.e. that requires the highest 
derivative on /) in fact comes from bounding (/, A s uj). 
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We summarise our results as: 
Theorem 2. Let the initial data v(0) G L 2 (M) and the forcing be bounded as 

(3.43) K s (f) := [V s+2 /fl + [V s d t /] < oo. 

Then there exist T s (\v(0)\ L2 , K s ; /x) and M S (K S ; /x) smc/i £/ia£ 
|V s w(i)| 2 2 QM S , 

(3.44) , /•«+*' 

^e-^t+O J e VT \V s u{r)\ 2 At <eM s 

for all t > T s . 

3.3. Higher-order Bounds. As in [13], one can obtain bounds that scale as e"/ 2 
for Cj when the force / is independent of time; see [1] . 

4. Stability and the Global Attractor 

When the forcing / is independent of time, the existence of the global attractor A 
follows, just as for the non-rotating 2d Navier-Stokes equations, from the uniform 
long-time bounds in Section [2] where the planetary rotation does not appear at all. 
In the non-rotating case, the Hausdorff dimension of A is bounded by 

(4.1) dim H ^<cG 2/3 (l+logG) 1/3 
where in our notation the Grashof number is 

(4.2) GiHV-VIWm 2 - 

The rotation not posing any extra essential difficulty, the usual analysis, e.g. [5j 
§9.2], carries over essentially line-by-line to our case, giving the bound (|4.ip also 
for the rotating case (|1.7|) . 

As discussed in the introduction, and following our results that the flow becomes 
more zonal ("ordered") as e — > 0, we expect the dimension of the attractor to 
decrease as e — > 0. In this section, we use a simple computation similar to that 
used for Theorem Q] to show that dim H A — for e sufficiently small. 

Theorem 3. Let the forcing f be time independent, dtf — 0, and assume the 
hypotheses of Theorem{j\ i.e. v(0) G L 2 (M) and 

(4.3) |V 2 /| i2 < ex). 

Then there exists an e»(|V 2 /|; /i) such that, for all e < £*, 

(4.4) dim H A = 0. 

Since A is connected, (|4.4[) implies that A consists of a single point, that is, 
a steady flow to which all bounded solutions converge. Following Theorems Q] 
and[2l this steady flow is nearly, but not exactly, zonal (except in the non-generic 
case when / = 0). Heuristically, an approximation to is the steady flow 

(4.5) = -»- 1 A- 1 f + eL- 1 f, 
which satisfies 

(4.6) -Lu^ +B(ul 1) ,u)i 1) )+iiAwi 1) = f 

up to 0(e). More careful work would be needed to determine u;* exactly. 
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In turbulence parlance, the smallness of e demanded by Theorem [3] implies that 
the Rhines scale [14] is so large that it overwhelms the entire spectral range, ren- 
dering the dynamics trivial. 

A general result related to ours is described in [3j ch. 18], where the trajectory 
attractor A e of a dynamical system depending on t/e (formally, in our case A e would 
simply be the attractor A for e > 0) converges weakly to the attractor Ao of the 
corresponding averaged system. Formally averaging our equations following this 
construction (which does not apply directly to our case, in which the oscillations 
have an infinite number of frequencies which accumulate at zero), we obtain purely 
zonal NSE, whose dynamics is trivial and whose attractor thus has dimension zero. 
This is of course consistent with our results: strong convergence at finite e of A to 
a point (which becomes zonal as e — > 0). 

PROOF. Fix a solution uj{t) of (|1.8|) that lives on A, so the bounds (|3.44j) hold for 
all t. We consider a nearby solution tu{t) + <f>(t). The linearised evolution equation 
for <f> is then 

d t <f> = -(V^A^w) • \7<f> - (V^A" V) ■ Vw(i) - -^A -1 ^ + ^A0 

(4.7) 

= -B(uj, 6) - B((j), uj) - -L(j> - fiA(f> =: C(tU. 

£ 

Multiplying this by <f> in L 2 and noting that (B(lu, <£>), <fi) = 0, we obtain 

(4.8) ^ l</,|2+/i|W|2 = 

= (fl(^),w) + (B(0,0),w). 

For the first term, we split <fi = ip + <f> in analogy with ui = Q + Co to get 

(4.9) (B(60),w) = (B(^),w) 

using the (now familiar) facts that B((f>, <p) — and all tilde-bar-bar terms vanish. 
Using Poincare inequality in (|4.8|) gives us [cf. (|3.10[) ] 

(4.10) |(e"^| 2 ) 
which integrates to 



(4.10) ^(e^|0| 2 )+ M e^|V0| 2 <2e^(i?(0,0),w) + 2e^(i?(^,^),^, 



\<j>{t)\ 2 e vt + v f |V0| 2 e^ r dr 
(4.11) Jo 

<|^(0)| 2 + 2 [ {(B(l^,Q) + (B^,4>),u>)}' 



o 

We bound the last term of the integrand using 
(B^,<P),Q) < c|V-V|l~|V0Uw| L 2 

(4.12) 

<c 4 |V0| 2 |^| L2 . 
The other term needs to be integrated by parts, 



(4.13) 



(B(<t>, j>), u) e UT dr = e (B n (0, $), u)(t) e ut - e (B n (4>, $),u)(Q) 

{v {B n (4>, $),w) + (B n (0, 0), d T co) + 2 (B n (6?0, 0), «)} e" T dr 
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(4.15) 
(4.16) 



where d^cf) = —B(uj, <p) — B((j>,w) — fj,A<j>. We bound the endpoint terms using 

(4.14) 2\(Bn(M),u)\ < c 5 |0| 2 |w'| ioo . 

It remains to bound the integrand in (|4.13|l : 

\{B n (4>,4>),u>)\ < c\d y 4>\ L2 \4>\ Li \oj\ L 4 

<c|V0| 2 |w| L4 
\(B n (4>,4>),d t u>)\ < c\d y 4>\ L 2\4>\ L i \d t u\ LS/ 2 

<c\V4>\ 2 \dM L , /2 

Recalling (j4.7|) for the last term in (|4.13|) . we bound 

\(B a {B(ct>,u>),4>),Q)\ < c|V-Vli-|Vw| i3 |^| il o|u / | i , /a 
<c|V0|V| L5/2 |Vw| i2 

\(B n (B(Lj,(t>),$),Q)\ < c\V- 1 u\ Lee \V^\^\ L io\u , \ L B / , 
(4.17) <c|V0|V| L5/2 |Vw| i2 

|(B n (A<^),u>)| = |(B n (V&V(?),w)| + \(B u (d y 4>,4>),u>')\ 

<c\Vt\ 2 \u'\ L „+c\Vt\ L S\ L M"X^ 
<c\V4>\ 2 \u>"\ L5/2 . 
Collecting, (|4.11|) now implies 

|«^)| 2 (l - ecsfi'®]^) + -sN(r) - c 4 \Q(r)\ L2 } |V0| 2 e"^ dr 

<e-^|^0)| 2 (l + ec 5 |^(0)| Lo .) 

where 

(4.19) N(t) := c 6 {fi\uj"\ L5/2 + |£DV/»|Vw| ia + |ftw| 18/a + |c3|i«}(t). 

By p. lip . / £ i/ 2 implies that u> £ H 3 with uniform bound in t since we are 
already on the attractor, and by Theorem [1] we can find an e* so small that, for 

£ < £*, 

(4.20) sup{dV(£) + c 4 |w(*)| L2 } < /i. 

i>0 

Requiring furthermore that e* also implies, for all e < £*, 

(4.21) ec 5 sup |w'(t)|i« < 1. 

t>o 

These and (|4.18l) then imply that 

(4.22) |0(t)| 2 < C(- • - )e-^|0(O)| 2 , 

in other words, all phase space volumes contract and thus the global attractor has 
dimension zero. 

It is clear from the above proof that our solution w(t) is linearly stable. Since 
(|4.7p only differs by B(<f>, <fi) from the nonlinear system, the fact that (B(<p, 0), <fi) — 
implies that stability also holds under the same hypotheses for the full nonlinear 
system. 
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Appendix A. An L°° Inequality 
Lemma 1. Let u and v e _ff 2 (M) have zero integrals and are L 2 orthogonal, 
(A.l) (u,v) L2 =0, 

and let W = U + v. Then the following Agmon inequality holds, 

/ I Awl \!/2 

(A.2) \u\ LOO <c\Vw\ (log-^P- + V 



c |Vw| 

Before the proof, we note that the interpolation inequality 

(A.3) |Vw| 2 < c 9 \w\ \Aw\ 

can be written as 

2 log |Vu>| < log \w\ + log |Aw| + logc 9 

^ A4 ^ & log|Vw| - log|w| < log|Aw| - log|Vw| +logc 9 
c \w\ c \Vw\ 

which can be used to simplify, e.g., |u)| ioo |Vw| i0 o when bounded using (IA.2I) . 

Proof. For most of this proof, up to (|A.7|) below, we follow [5J Lemma 7.1] 
exactly. For conciseness, we put L\ = Li = 1 but keep the Poincare constant Co- 
With k > 0, we expand u in Fourier series 

(A.5) u{x) = E\ k \< K u ke ik X + E| fc |>« u k e ik x =: u<{x) + u>{x), 
and analogously for v and w. Then 

\u\ Lac = SU Pa . \j2k u k^ k ' X \ < E|fc|< K l«fcl +E|fc|> K l«J 

=--E < |fer 1 |fc|KI + E > |fer 2 |fe| 2 KI 
(A.6) w^<i,.,-aV>V^< .«..*... , 2 V /2 



^(E^r 2 ) (e^^k 



f \ 1/2 / \ 1/2 

E^fcr 4 fe > i*n«* 19 



Now on the right-hand side, E < 1^1 2 — c l°g' K and E > 1^1 4 — c / r2 7 so fixing 
(A.7) k = |A«;|/(co|V«;|), 
the lemma follows from 



(A.8) V Co|Vu;|y 

|Vw|(log- 



I Aw I \ 1 /2 
< clVwII log ' — ' , +c|Vw 
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